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ABSTRACT 



£f>l . We have undertaken an HST Space Telescope Imaging Spectrograph survey of 54 late type spiral galaxies to study the scaling relations 

between black holes and their host spheroids at the low mass end. Our aim is to measure black hole masses or to set upper limits for a 
sizeable sample of spiral galaxies. In this paper we present new Space Telescope Imaging Spectrograph (STIS) observations of three 
spiral galaxies, NGC 4303, NGC 3310 and NGC 4258. The bright optical emission lines Ha X 6564A, [Nil] XX 6549, 6585A and 
[SII] XX 6718, 6732A were used to study the kinematics of the ionized gas in the nuclear region of each galaxy with a ~ 0.07" spatial 
resolution. Our STIS data for NGC 4258 were analyzed in conjunction with archival ones to compare the gas kinematical estimate of 
the black hole mass with the accurate value from Z^O-maser observations. 
In NGC 3310, the observed gas kinematics is well matched by a circularly rotating disk model but we are only able to set an upper 
limit to the BH mass which, taking into account the allowed disk inclinations, varies in the range 5.0 x 10 6 - 4.2 x 10 7 M o at the 95% 
' confidence level. 

In NGC 4303 the kinematical data require the presence of a BH with mass M BH = (5.0)^ ^ X 10 6 M o (for a disk inclination i = 70 
deg) but the weak agreement between data and disk model does not allow us to consider this measurement completely reliable. If the 
allowed inclination values are taken into account, M BH varies in the range 6.0 x 10 s - 1.6 x 1O 7 M at the 95% confidence level. 
In NGC 4258, the observed kinematics require the presence of a black hole with M BH = (7.9)+^ 2 x 10 7 M o (i = 60 deg) and, taking 
into account reasonable limits for the inclination, M BH is in the range 2.5 x 10 7 - 2.6 x 10 8 M o at the 95% confidence level. This result 
is in good agreement with the published value (3.9 ± 0.1) x 10 7 M G , derived from H^O-maser observations. As in the case of NGC 
4303, the agreement between observed and model kinematics is not strong but this does not affect the recovery of the correct M BH 
value. 

Our attempt at measuring BH masses in these 3 late type Sbc spiral galaxies has shown that these measurements are very challenging 
and at the limit of the highest spatial resolution currently available. Nonetheless our estimates are in good agreement with the scaling 
relations between black holes and their host spheroids suggesting that (i) they are reliable and (ii) black holes in spiral galaxies 
follows the same scaling relations as those in more massive early-type galaxies. A crucial test for the gas kinematical method, the 
correct recovery of the known BH mass in NGC 4258, has been successful. 

Key words. Black hole physics - galaxies: spiral - galaxies:bulges - galaxies:kinematics and dynamics - galaxies:nuclei - galax- 
ies:individual:NGC 3310- galaxies:individual:NGC 4303- galaxies:individual:NGC 4258 



1. Introduction 

It is becoming increasingly clear that the energy output from 
supermassive black holes (BH) at galaxy centers plays a pivotal 
role in the formation and evolution of galaxies. 

Send offprint requests to: Guia Pastorini 

* Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. These observations are asso- 
ciated with proposal 8228. 



It is now widely accepted that the mechanism powering 
Activ e Galactic Nu clei (AGN) is accretion of matter onto a 
BH (ISarpetenTl 964). with masses in the 10 6 - 10 10 M© range 
dFerrarese & Fordl2005l) . This paradigm, combined with the cos- 
mological evolution of AGNs, implies that most, if not all, local 
galax ies should host a BH in their nuclei (e.g. iMarconi et al.l 
2004 and references therein). 

B H are detected in about 40 galaxies, mostly early-types 
(E/SO: IFerrarese & Fordll2005l and references therein). The BH 
mass (Mbh) is closely related to the stellar velocity dispersion 
(o-.: IFerrarese & Merrittll20rjol [Gebhardt etai] |200of) and to the 
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host bulge mass or luminos ity (L sp h, iKormendv & Richstond 
119951: IMarconi & Hunlll2003"l) . These correlations reveal a tight 
link between the growth of central BH and the formation and 
evolution of the host galaxy. These lo cally detected BHs were 
probably grown during AGN activity (Mar coni et al.|[2004h and 
during these phases they had a major impact on the host galaxy: 
the feedback from the accreting BH (i.e. from the AGN) can be 
so strong as to suppress star formation in the host galaxy and 
is possibly the physical mechanism responsible for setting the 
close relat ions between BH and host galaxy structural parame- 
ters (e .g. lMenci et alj2003tlGranato et alj2004tlDi Matteo et alj 
120051 1 

The number of BH detections in spiral galaxies is small (cur- 
rently 7, iFerrarese & Fordll2005l) and it is clearly mandatory to 
increase this number in order to understand how common BHs 
are in spirals and if they follow the same correlations as early 
type galaxies. Even constraining upper limits on Mbh could be 
considered as a significant improvement given the scarcity of 
measurements and the difficulties in performing the required ob- 
servations in spiral galaxies. Moreover, M^n-galaxy correlations 
are studied only in the 3 x 10 6 - 3 x 1O 9 M mass range and it 
is expected that spiral galaxies might host low mass black holes 
below 10 7 M G . 

The direct methods adopted to measure Mbh in the nearby 
universe use gas or stellar kinematics to gather information 
on the gravitational potential in the nuclear region of the 
galaxy. The stellar kinematical method has the advantage that 
stars are present in all galactic nuclei and their motion is al- 
ways gravitational. The drawback is that it requires relatively 
long observation times in order to obtain high quality ob- 
servations and that stellar dynamical mod els are very com- 
plex, potentially plagued by indeterminacy dValluri et alj r2004: 
ICretton & Emselleml |2004|) . Conversely, the gas kinematical 
method is relatively simple, it requires relatively short observa- 
tion times for the brightest emission line nuclei, even if not all 
galaxy nuclei present detectable emission lines. Another impor- 
tant drawback is that non-circular or non-gravitational motions 
can completely invalidate this method. Since the observed corre- 
lations are based on BH masses obtained with different methods, 
it is important to check whether these methods provide consis- 
tent and robust results. A comparison between Mbh measure- 
ments obtained with gas and stellar kinematics so far has been 
perfo rmed only for the Milky Way dGenzel et al.ll2000h iGenzeTI 
120061) and Centaurus A ( IMarconi et alj|2005t ISilge etal.|[2005h . 
where the two independent estimates are in excellent agreement. 
This comparison has b een attempted al s o for IC 1459, NGC 
4335 and NGC 33 79 dCappellaril l200l IVerdoes Kleiin et al. I 
I2002t IShaoiro et alJ^OOaT but in those cases the gas kinemat- 
ical method does not provide reliable MBH estimates because 
the observed velocity field is not reproduced by a circularly ro- 
tating disk. 

All methods used to directly detect BHs and measure their 
masses require high spatial resolution in order to resolve the BH 
sphere of influence, i.e. the region where the gravitational in- 
fluence of the BH dominates over that of the host galaxy. The 
radius of the BH sphere of influence is traditionally defined as 
rBH = GMbh/v 2 , where G is the gravitational constant and cr is 
the stellar velocity dispersion, and is usually < 1". Very high 
spatial resolution is provided by the Hubble Space Telescope 
(HST) whose angular resolution is approximately 56 = 0.07" 
(Full Width at Half Maximum, FWHM, of the Point Spread 
Function, PSF) at A ~ 6500A. Let us assume that a spiral galaxy, 
located at a distance of 20 Mpc, hosts a BH with M B h - 10 7 M p 
and follows the Mbh ~ cr correlation (Tremai ne et alj 12002): 



the expected value of cr is then ~ 105 km/s. The radius of the 
black hole sphere of influence expected for this object is thus 
rBH ~ 0.04" very close to the HST angular resolution, demon- 
strating how challenging these observations are. 

We have undertaken a spectroscopic survey of 54 spirals us- 
ing the Space Telescope Imaging Spectrograph (STIS) on HST. 
Our sample was extracted from a comprehensive ground-based 
study by D.J. Axon et al. (in preparation) who obtained Ha 
and [Nil] rotation curves, at seeing-limited resolution of 1", of 
128 Sb, SBb, Sc and SBc spiral galaxies selected from RC3. 
By restricting ourselves to galaxies with recession velocities 
V < 2000 km s _1 , we obtained a volume-limited sample of 
54 spirals known to have nuclear gas disks. The sample spans 
a wide range in bulge mass and concentration. The systemic ve- 
locity cutoff was chosen so that we can probe close to the nuclei 
of these galaxies and detect low-mass BHs. 

This paper is par t of a series reportin g the results of our HST 
program. In Paper I ( Ma rconi et al.ll2003l) we reported the mod- 
eling of the central mass concentration in NGC 4041, the first 
galaxy to be observed as part of the program. Paper I also con- 
tained a detailed description of the modelin g techniques that we 
use to determine BH masses. In Paper II faughes et alJl2003h 
we presented the atlas of STIS spectra obtained from the suc- 
cessful observations and we created color maps from the optical 
STIS images when archival Near-Infrared Camera and Multi- 
Obje ct Spectrometer (N ICMOS) images were available. In Paper 
III dScarlata et a l. 2004) we presented the STIS images acquired 
in the program and analyzed surface brightness p rofiles derived 
from the images. In Paper IV ( iHughes et al]|2005h we used both 
color information from these images and the spectra presented 
here to in vestigate the age of t he central stellar populations. In 
Paper V (Atk inson et alj|2005l) we reported the BH mass mea- 
surements in two Sbc galaxies of the sample, NGC 1300 and 
NGC 2748. Here (Paper VI) we present the analysis of three ad- 
ditional Sbc galaxies from our sample: NGC 3310, NGC 4303 
and NGC 4258, selected from the remaining sources because 
they combine high S/N spectra with evidence for rotation. 

NGC 3310 (Arp217, UGC 5786), located in the vicinity of 
the Ursa Major clu ster, is an Sbc galaxy ("S AB(r)bc pec" in the 
RC3 catalogue bv lde Vacouleurs et al.ll 199 lb with an inclination 
of the galactic disk of about i ~ 40 deg ( Sanchez-Porta feTaTI 
2000). It is a very well known starburst galaxy char acterized 
by a disturbed morphology (see Elme green et alj|2002l for a re- 
cent review) which suggests that star formation has been trig- 
gered by a collision with a dwarf companion during the las t 
~ 10 8 yr dKregel & Sancisil l200ll: iBahck & HeckmSl [l98l 
Its central starbust region dDiaz et alj 120001 : lElmegreen et alj 
2002; Chandar et al. 2005), as well as its young stellar clusters 
dde Griis et alj 2003a b). have been studied in detail revealing 
the presence of very young stars (from ~ 10 7 to 10 8 yr). The 
bright inner region is dominated by a two-armed open spiral pat- 
tern observed in Ha and the inner part of this well-developed 
pattern connects to a ~ 900 pc dia meter starburst ring, sur- 
round i ng the blue compact nucl eus dvan der Kruit & de Bruvnl 
119761 : iBalick & Heckmanl 1 198 lh . The gas kinematics in the 
galactic disk is disturbed, characterized by non-circular mo- 
tions and strong streaming along the spiral arms, strengthen- 
ing the idea of a recent merger event (Rregel & Sancisil l2001; 
Muld er & van Drielll996l) . Concerning the nuclear region, the 
rotation center of the gas is offset with respect to the stellar 
continuum isophote s by ~ 1.5" ± 0.3" toward PA 142 deg 
dvan der Kruitl 1976b . The STIS acquisition images of NGC 3310 
obtained with our HST project indicate th e presence of a stron g 
dust lane crossing the center of the galaxy dScarlata et alj |2004). 
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From the Lyon/Meudon Extragalactic Database^ (LEDA) NGC 
3310 has an average heliocentric velocity of 1216 km/s, after 
correction for Local Group infall onto Virgo. Assuming Ho = 70 
km s" 1 Mpc -1 , this implies a distance of D ~ 17.4 Mpc (i.e. 
l"~84pc). 

The late-type Sbc galaxy NG C 4303 (M61) is a d ouble- 
barred galaxy in the Virgo Cluster (Binggelli et aTll 19851) at dis - 
tance of 16.1 Mpc (1" ~ 78 pc, see e.g. lSchinnereret al. 2002). 
A large-scal e bar of about 2.5 kpc (~ 35") lies inside the outer 
spiral arms (lLaine et al.ll2002t lErwinll2004 . In addition to the 
large-scale bar, this galaxy hosts a second inner bar of about 0.2 
kpc (~ 2") length that is surrounded by a c ircumnuclear starburst 
ring/spiral of ~ 0.5 kpc (~ 6") diameter dColina & Wadall2000l 
Perez-Ram irez et al.ll2000t IColina et al.lll997l) . Based on the op- 
tical emission li ne ratios it is classified as a Seyfert 2/L INER 
borderline AGN (Filippenko & Sargent 1985; Colina & Arribas 



1999) and, in fact, it constitutes a nice example of the starburst- 
AGN connection. A young (~ 4 Myr), massive stellar cluster and 
a low-luminosity active nucleus seem to coexist within the cen- 
tral 3 pc (~ 0.04") of the galaxy with the starburst contribution 
being the dominant one in the UV spectral region dColina et alj 
2002). The AGN pres ence is revea led through the luminous near- 
infrared p oint source dColina et al.l2002l) and through hard X-ray 
emission ( Jimenez-Bail on et al .120031) . The gas kinematics of the 
galactic disk is affected by the presence of the large bar but, in 
the nuclear region, the gas seems to show ordered ro tation. From 
integral field observations. lColina & Arribasld 19991) find that the 
H/3 velocity field at ~ 1.5" resolution is consistent with a circu- 
larly rotating massive disk of ~ 300 pc radius (~ 4"). CO obser- 
vations confirm this behavior and indicate that, at a few arcsec 
resolution, the gas within the central 0.7 kpc (~ 10") is follow- 
ing almost pure circular rotation with the dyna mical center lo- 
cated within 1" from the ne ar-IR point source (Koda& Sofue 
2006; Schinnerer et al. II2002I) . At 2" resolution. lSchinnerer et alj 
(2002) find the presence of an asymmetric spiral pattern inside 
a radius of ~ 400 pc (~ 5") which is consistent either with a 
self-gravitating m = 1 mode (hence gas kinematics would be 
unaffected by the small inner bar) or with a ring of star forma- 
tion heavily extincted north of the nucleus. 

NGC 4258 (Ml 06 ) is a SABbc nearby gal axy (D = 7.2 + 0.3 
Mpc, i.e. 1" ~ 35 pc iHerrnstein et all 19991) spe ctroscopicallv 
classified as a 1.9 Sevfert galaxy fao et al.lll997l) . This galaxy 
is very well known because it represents the b est case for a 
nuclea r supermassive BH after the M ilky Way . Mivoshi et alj 
d 1995b . iGreenhill et all d!995allbl) and IHerrnstein et al.1 d 1999b 
studied the H2O maser line emission of a rotating gas disk near 
the galaxy center. This disk is nearly edge-on (i — 81 ± 1 deg) 
and its kinematics is nicely explained by circular rotation around 
a (3.9+0. l)xlO 7 M mass confined within 0.1 3 pc of the nucleus, 
indicating very compelling evidence for a supermassive nuclear 
black hole. The active nucleus has bee n detecte d through radio 
dTurner & Holll994t ICecil et al.|[l995l). X-ray dFruscione et al.l 
120051; IWilson et al.M2001t iMakishima et al.lll994l) and in frared 
observations dCharv & Becklinl 119971; ICharv et al.l 120001) . The 
strong polarization of the relatively broad optical emission 
lines further suppor ts the existence of an obscured active nu- 
cleus in NGC 4258 dWilkes et al.lll995l) . In particular, the AGN 
is characterized by a rad io jet which pro pagates perpendic- 
ularly to the maser disk dCecil et al.l 11995b . On large scales, 
"anomalous arms" were di scovered through Ha imaging by 
ICourtes & Cruveliierl dl96ll) . which noted their unusual diffuse 
and amorphous appearance, in strong contrast with the knotty 
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structure o f normal spiral arm s due to the presence of HII- 
regions (se elWilson et al.l2001lf or a recent review on the anoma- 
lous arms). IWilson et all (|2001|) interpreted the anomalous arms 
as jet-cloud interaction between the radio jet and the galactic 
disk: the jets shock the normal interstellar gas along the first 350 
pc (~ 10") of their length, causing the hot, X-ray emitting co- 
coons surrounding the radio jet. To summarize, in the nuclear 
region of NGC 4258, the radio jet propagates perpendicularly to 

the maser disk at PA 5 deg and the shock edge between the 

hot X-ray cocoon and the dense gas in the galactic disk is at PA 
-33 deg. 

The aim of this paper is to present the gas kinematical analy- 
sis of these three spiral galaxies in order to measure the mass of 
their nuclear BHs. The case of NGC 4258 is particularly signifi- 
cant because this galaxy is an important test for the gas kinemat- 
ical method since its BH mass is accurately known from high 
spatial resolution maser observations described above. 

The structure of the paper is as follows. In Sec.|2]we present 
our HST/STIS and archival HST/NICMOS-ACS observations as 
well as the data reduction procedures. In Sec. |3.1| we describe the 
kinematical analysis and we present the kinematical parameters 
measured from our data. In Sec. l3.2l we describe the method fol- 
lowed to obtain the stellar contribution to the gravitational po- 
tential for each galaxy. In Sec. [4] we describe the kinematical 
modeling and how it is possible to measure Mbh- In Sec. [5] we 
present the results for all galaxies. These results are discussed 
and compared with the expectations from the Mg^-galaxy cor- 
relations in Sec. [6] and, in Sec. [7] we summarize our work and 
present our conclusions. 



2. Observations and data reduction 

For each galaxy we analyze our own HST-STI S spectra and 
archiv al HST-NICMOS images (Fig. [Q see also iHughes et all 
2003). In the case of NGC 3310 we consider also an archival 
narrow band Ha image acquired with HST-ACS (Fig. [2j. In the 
case of NGC 4258 we combine our HST-STIS data with archival 
ones. In Tab.[TJwe present a list of all available observations and 
the following subsections are dedicated to an accurate descrip- 
tion of the analysis of spectra and images. 

2.1. STIS Observations 

The three galaxies were observed with STIS onboard the HST 
The adopted observational strategy consisted in obtaining spec- 
tra at three parallel positions with the central slit centered on 
the galaxy nucleus, and the flanking ones offset by 0.2". The 
archival STIS observations of NGC4258 were obtained with a 
similar strategy but the offset of the flanking slits was 0.15". 
These different slit positions are shown in Fig. [TJ and labeled 
'NUC, 'OFF1' and 'OFF2'. The log of the observations is pre- 
sented in Tab. Q] 

At each slit position, we obtained 2 spectra with the galaxy 
nucleus placed in different positions along the slit, shifted by an 
integer number of pixels, in order to remove cosmic-ray hits and 
hot pixels. Each spectrum was obtained with the G750M grating 
and the 0.2" slit which provided a dispersion of 6 A = 1.1 08 A 
pix -1 and a spectral resolution of R = A/(26A) ~ 3000. In the 
case of the archival NUC position in NGC4258, the slit size was 
0.1". All the spectra, except for the NGC4258 archival and our 
nuclear ones, were obtained with a 2x2 on-chip binning resulting 
in a pixel spatial scale of 0.101" pix -1 . See Tab. Q] for details on 
slit size and binning. 
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Fig. 1. Left panel: HST-NICMOS F160W image of NGC 3310 with the three different parallel positions of the slits used for spectra 
acquisition. North is up and East is left. The gray-scale was chosen to have darker points for higher intensity values. The slit 
positions used in the observations are labeled as NUC, OFF1 and OFF2 (see text for more details). Central and right panels: NGC 
4303 and NGC 4258, notations as in left panel. 
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Fig. 2. HST-ACS narrow band Ha image of NGC 3310 with slit positions overlaid. North is up and East is left. 



The data were first redu ced using the STIS pipeline, cal- 
stis dKim-O uiiano et al. 2003), and the most up-to-date reference 
files. The raw spectra were first reprocessed through the calstis 
pipeline and standard pipeline tasks were used to obtain flat-field 
corrected images. The two exposures taken at a given slit posi- 
tion were then realigned with a shift along the slit direction (by 
an integer number of pixels), and the pipeline task ocrreject was 
used to reject cosmic rays and hot pixels. Subsequent calibration 
procedures followed the standard pipeline reduction describ ed in 
the STIS Instrument Handbook dKim-Ouiiano et ai] |2003): i.e., 
the spectra were wavelength calibrated and corrected for two- 
dimensional distortions. The expected accuracy of the wave- 
length calibration is 0.1-0.3 pixels within a single exposure and 
0.2-0 .5 pixels among different exposures (Ki m-Ouiiano et all 
120031) . which converts into~ 3-8 km s _1 (relative) and ~ 5 - 13 
km s" 1 (absolute). The relative error on the wavelength calibra- 
tion is negligible for the data presented here because our analysis 
is restricted to the small detector region including Ha, [Nil] and 
[SII] (AA < 200 A). 



2.2. Images 

HST- NICMOS images obtained with Camera 2 dMacKentv et alj 
\\99% and the F160W filter (~ H band) were retrieved from the 
HST archive. Their spatial scale is 0.076"/pix and they have a 
field-of-view of ~ 20 X 20". The images, with the slit positions 
superimposed, are shown in Fig. Q] and their det ails in Ta b. [2] 
Imagin g data reduction is described in detail by Hughes et al. 
(120031) . 

For NGC 3310 we have also retrieved a narrow- band filter 
(F658 N) image obtained with ACS (camera WFC1) dMack et al.l 
2003 ). The imag e was recalibrated using the ACS pipeline 
dMack et alj2003l) with the most up-to-date calibration files. The 
final spatial scale of the calibrated image is 0.05"/pix and its to- 
tal field of view is ~ 300" x 300". The image is displayed in 
Fig. [2] and its details in Tab. [2] 



3. Results 

In this section we present the results from the spectroscopic and 
imaging observations and the preliminary steps required to per- 
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ID 


PID 


Rootname 


Date 


9 [degf 


t exp [s] 


Offset ["] 


Slit size [" | 


Binning 


NGC 3310 


NUC 


8228 


o5hl7030,40 


Feb 11, 2000 


170 


864 


0.0 


0.197 


2x2 


NGC 3310 


OFF1 


8228 


o5hl7010,20 


Feb 11, 2000 


170 


480 


-0.2 


0.197 


2x2 


NGC 3310 


OFF2 


8228 


o5hl7050,60 


Feb 11, 2000 


170 


576 


0.2 


0.197 


2x2 


NGC 4303 


NUC 


8228 


o5h735030,40 


March 31, 2000 


135 


720 


0.0 


0.197 


2x2 


NGC 4303 


OFF1 


8228 


o5h735010,20 


March 31, 2000 


135 


480 


-0.2 


0.197 


2x2 


NGC 4303 


OFF2 


8228 


o5h735050,60 


March 31, 2000 


135 


576 


0.2 


0.197 


2x2 


NGC 4258 


NUC 


8228 


o5h734030,40 


August 01, 2000 


246 


864 


0.0 


0.197 


lxl 


NGC 4528 


OFF1 


8228 


o5h734010,20 


August 01, 2000 


246 


576 


-0.2 


0.197 


2x2 


NGC 4258 


OFF2 


8228 


o5h734050,60 


August 01, 2000 


246 


576 


0.2 


0.197 


2x2 


NGC 4258 


NUC 


8591 


o67104010,20,30 


March 16, 2001 


146 


144 


0.0 


0.095 


lxl 


NGC 4528 


OFF1 


8591 


o67 104040 


March 16, 2001 


146 


158 


0.15 


0.197 


lxl 


NGC 4258 


OFF2 


8591 


o67 104050 


March 16, 2001 


146 


120 


-0.15 


0.197 


lxl 



8 Position angle of the slit (N to E). 



Table 2. Log of imaging observations 





PID 


Rootname 


Date 


Instrument 


Camera 


Filter 


6(°y 


texpis) 


NGC 4303 


7330 


n3zl32u010 


April 03, 1998 


HST-NICMOS 


NIC2 


F160W 


-103.42 


320 


NGC 3310 


7268 


n49j06tkq 


November 12, 1997 


HST-NICMOS 


NIC2 


F160W 


80.53 


128 


NGC 3310 


9892 


j8og20010 


October 21, 2003 


HST-ACS 


WFC1 


F658N 


0.014 


950 


NGC 4258 


7230 


n46801040 


November 21, 1997 


HST-NICMOS 


NIC2 


F160W 


85.54 


56 



8 Position angle of image y axes (deg e of n). 



form the kinematical modeling which can result in a BH mass 
estimate. In particular in Sec. 13.11 we first describe the proce- 
dure adopted to estimate the kinematical moments along the slits 
(surface brightness, velocity, velocity dispersion and asymmetry 
parameters) and we then present the results for each galaxy. In 
Sec. l3.2l we briefly describe the procedure used to derive the con- 
tribution of the stellar mass to the gravitational potential and we 
then present the results for each galaxy. 

3.1. Kinematical analysis 

The emission lines analyzed to extract the kinematical pa- 
rameters are Ho- i6564A, [Nil] /U6548, 6583A and [SII] AA 
6716,6731A. Line-of-sight velocities, velocity dispersions, and 
surface brightnesses along each slit were obtained by fitting the 
emission lines in each row of the continuum-subtracted two- 
dimensional spectra. In all cases the emission lines were mod- 
eled with a Gaus s-Hermite series ( Ivan der M arel & Franx 1993; 
Cappellari 2004): 



4 

1 + 2j h mH m (y) 

m=3 



(1) 



where v is the velocity, y = (v - vo)/<x and H m with m = 3, 4 
are the Hermite polynomials of order 3 and 4, respectively 
(Cappellari 2004). The best fit parameters (line flux, average ve- 
locity, velocity dispersion, ^3 and h\) we re determined by min- 
imizing a penalized x 2 as described in Cappellari 2004] The 
Gauss-Hermite expansion is usually adopt ed in the ana l ysis of 
stellar kinematics from absorption lines; iMarconi et alj ([2005) 
provide a complete description of the application of this method 
to the analysis of emission line spectra. 

During the fit the only constraints we imposed were to link 
the kinematical parameters of lines belonging to the same dou- 
blet. Lines from the [Nil] or [SII] doublet were constrained to 
have the same velocity, velocity dispersion and /13, /14. In the 
case of [Nil] we also imposed a fixed flux ratio of 1 : 3 between 



the faintest and the brightest line of the doublet. The velocity, 
velocity dispersion, I13 and /14 for Ha and for the [Nil] and [SII] 
doublets were then left free to vary during the fit. Finally, for 
each slit position we considered the flux weighted average of all 
kinematical quantities. 

Fig.[3]and Fig.|4]present the resulting kinematical quantities 
measured along all slit positions for NGC 3310 and NGC 4303. 
The kinematical analysis of NGC 4258 followed a slightly dif- 
ferent procedure and is presented at the end of this section. 

The use of the Gauss-Hermite expansion allowed us to repro- 
duce asymmetries in the emission line profiles in a more robust 
way than, e.g., the use of t wo Gaussian functions (see the dis- 
cussion in lCappellarill2004l) . In the case of NGC 3310 and NGC 
4303 the coefficients of the Hermite polynomials /13, /14 are small 
(< 0.1 - 0.2) so that the fitted emission line profiles are similar 
to pure Gaussian functions and, indeed, the rotation curves ob- 
tained by imposing /13 = h\ — are not significantly different, 
only less smooth. 

Both NGC 3310 and NGC 4303 are characterized by very 
peaked surface brightness distributions which suggest the exis- 
tence of two different regions: a nuclear one, identified by the 
peak in surface brightness, and an extended one characterized by 
a surface brightness which is much smoother and more diffuse. A 
similar situation was found in NGC 4041 (Ma rconi et afll2003l) 
and is consistent with the findings bv lKnapenl d2005l) at ground- 
based spatial resolution. In Fig. [Hand Fig.[3]the 'nuclear' region, 
identified on the basis of the surface brightness distribution (see 
Sec.|5]l, is delimited by vertical dashed lines. 

The velocity curves are different for the two galaxies. In 
NGC 3310, the nuclear rotation c urve has the typical S feature 
expected from a rotating disk (e.g. lMarconi et al. 2005) with an 
amplitude of ~ 200 km/s. On the contrary, the velocity curve 
of NGC 4303 seems less regular and the S shape appears more 
distorted. The amplitude is also much smaller, ~ 50 km/s. Both 
galaxies show a large scale velocity gradient, such as that ex- 
pected from the rotation of the galactic disk, but the rotation 
curves are disturbed. The velocity dispersions of the gas do not 
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Fig. 3. Observed kinematical moments for NGC 3310 obtained with the fitting procedure described in the text. The first row shows 
the average emission line surface brightness (units of l0~ l6 ergs~ l cm~~ 2 arcsec~ 2 ). The second row shows the velocity along the line 
of sight (units of km/s) and the third one the velocity dispersion of the line (units of km/s). From left to the right, the three columns 
refer, respectively, to the NUC, OFF1 and OFF2 slit positions indicated in Fig.Q] The vertical errors-bars, when visible, represent the 
formal errors of the fit, while the horizontal bars are the sizes of the bins used to extracted spectrum (in general one detector pixel). 
The vertical deshed lines delimit the nuclear region, which we have identified on the basis of the surface brightness distributions in 
the upper panels and which will be analyzed in the kinematical model. 



show the behavior expected from a rotating disk (smooth plateau 
and peak at the position of the BH) but are more or less irregu- 
lar suggesting that turbulent or non circular motions might affect 
the kinematics. 

In the case of NGC 4258 we adopted a modified procedure 
with respect to the one used for NGC 3310 and NGC 4303. In 
this galaxy the Ho- and [Nil] lines are severely blended and, be- 
ing a Seyfert 1 .9, one has to take into account the possible pres- 
ence of a broad line. 

We first assume that Ha, [Nil] and [SII] share the same kine- 
matics. Although there might be differences due to the fact that 
Ha can be emitted also by gas at different stages of excitation, 
our assumption will in any case result in average kinematical 
quantities, weighted with the line surface brightness. We col- 
lapse the central 4 or 6 rows in order to obtain a high S/N nu- 
clear spectrum and we fit the [SII] doublet with the same con- 



straints adopted for NGC 3310 and NGC 4303. The kinematical 
parameters found from the [SII] doublet are then used to recon- 
struct the Ha and [Nil] narrow lines which are then subtracted 
from the spectrum. The flux of the Ha and [Nil] narrow lines 
is adjusted so that the residual spectrum does not show obvious 
structure due to under/over subtracted narrow lines. The resid- 
ual spectrum is that showing the broad line and any other com- 
ponents which have a different kinematics to that of the [SII] 
doublet (see Fig. |5). It is interesting to note that the broad Ha 
line (FWHM~ 1600 km/s) does not have a Gaussian profile but 
a pro file resembling that of double humped broad lin e galaxies 
(e. g. lEracleous & Halpern Il2003t IStrateva et al. 1120031 and ref- 
erences therein. ) [Nagao et all (120061) modelled the BLR of their 
quasar templates with a double power law function. Here, we 
find that a good model profile for the BLR is a three piece bro- 
ken powerlaw as shown in Fig. [6] 
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Fig. 4. Observed kinematical moments for NGC 4303. Notation as in Fig. [3] 




Fig. 5. The black line shows the observed NGC 4258 nuclear spectrum, the green one the narrow lines used for subtraction; the 
residuals are drawn in red. These residuals represent the Broad Ha Line. The left and the right panels refer, respectively, to archival 
and our own data. 



Note that the BLR profiles derived from our and archival 
spectra are consistent, within the uncertainties, in shape and nor- 
malization. If the broad line flux did not vary significantly from 
one observation to another, the above similarity implies that the 
NUC slits were well centered on the nucleus position. We then 
refitted all emission lines in nuclear spectra with the components 



(modeled with the Gauss-Hermite expansion and constrained to 
share the same kinematics) and a broad Ha modeled as described 
above. Except for the flux, all parameters of the broad Ha line 
derived from these fits were then kept fixed in the subsequent 
row by row fits. I13 and h\ were set to in the low S/N spectra of 
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Fig. 7. Example of line fitting of the Ha, [Nil] and [SII] lines in the nuclear region of NGC 4258. The left and right panels refer, 
respectively, to our and archival data. The lower window of each panel shows the residuals of the fit. 



NUC and in all spectra of OFFl,2. Fig. [7] shows an example of 
the fit for ours and the archival data. 

The kinematics extracted with this fitting procedure are 
shown in Fig. [8] In contrast with NGC 3310 and NGC 4303, 
the rotation curves are limited only to within 0.5" of the cen- 
tral source. As for the two other spiral galaxies, we detect the 
bright nuclear emission but we lack the fainter extended com- 
ponent. With the limited spatial coverage, the rotation curves 
appear more regular than in the cases of NGC 3310 and NGC 
4303. The velocity amplitudes of the rotation curves in our data 
are about ~ 200 km/s while these drop to ~ 100 km/s for the 
archival data. This is not surprising since our spectra and the 
archival ones are taken at different position angles. The veloc- 
ity dispersions, apart from some deviant points, seem to exhibit 
a nuclear peak at about ~ 200 - 250 km/s and much smaller 
constant values at the offnuclear positions ~ 150 km/s with a 
behavior similar to that of a rotating disk. 

A possible cause of concern for the accuracy of our kine- 
matical measurement is given by the fact that the point spread 
function of HST is undersampled by the pixels of the STIS CCD. 
When spectral images with unresolved components are rectified, 
interpolation along the undersampled spatial axis produces arti- 



facts in the rectified image. An example of the artifacts in the 
peak row of a G750M image is shown in F igure 11.9 (Section 
1 1.3.5 ) in the STIS Instrument Handbook dKim-Ouiiano et all 
2003). The spectrum in each row of the image undulates across 
the detector as flux is alternately retained, gained, or lost by 
the interpolation. Adjacent rows have undulations that are out 
of phase, since the interpolation inappropriately pushes flux out 
of one row and into the next. Undulations were obvious in the 
continuum spectrum of NGC4303. They presented no problems 
for the analysis of the emission lines, however. The lines were so 
narrow that the underlying continuum was well represented by 
a short linear segment. NGC 4258 was a somewhat more prob- 
lematic case, since it has a spatially unresolved spectral feature 
(a broad H a component) which is broad relative to the structure 
of the undulations. This undulating feature must be subtracted 
from the rows of the spectrum so that the parameters of the nar- 
row lines can be measured. The undulation structure cannot be 
directly observed in this case, since it is mixed in with the spa- 
tial and spectral structure of the H a + [Nil] blend. We assessed 
its impact by artificially dividing our nuclear spectra with sev- 
eral observed undulation patterns and remeasured kinematical 
parameters. The impact on the velocity measurements is in gen- 
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Fig. 8. Observed kinematical moments for NGC 4258: the first three columns refers to our data, the other ones represent the archival 
data. Archival and our own data were acquired with the slit placed at two different position angles as shown in FigQ] 



eral small. Moreover, the size of the residuals in the original fits 
sets limits to the magnitude of the neglected undulations. Our 
tests indicate that undulations of this magnitude would not sig- 
nificantly affect the velocity measurements. 

3.2. Light Profile 

In order to account for the contribution of the stars to the grav- 
itational potential in the nuclear region, the stellar luminosity 
density has to be derived from the observed surface brightness 
distribution. The inversion procedure to derive this density dis- 
tribution from the observed surface brightness is not unique if 
the gravitational potential does not have a spherical symmetry, 
such as can be revealed by circular isophotes. Assuming that 
the gravitational potential is an oblate spheroid, the inversion 
depends on the knowledge of the potential axial ratio, q, and 
the inclination of its principal plane with respect to the line 
of sight, i. As these two quantities are related by the observed 
isophote ellipticity, we are left with the freedom of assum- 
ing different galaxy inclinations to the line of sight. Following 
van der Marel & van den Bosch 1 (119981) . we assumed an oblate 
spheroid density distribution parameterized as: 



1 + 



(2) 



m is given by m 2 - x 2 + y 2 + z 2 /q 2 ; xyz is a reference sys- 
tem with the xy plane corresponding to the principal plane of 
the potential and q is the intrinsic axial ratio. When q — 1 the 
model is spherical and m corresponds to the radius. The above 
stellar density is integrated along the line of sight to provide 
the surface brightness distribution on the plane of the sky. This 
is then convolved with the ins trumental Point Spre ad Function 
(PSF) (derived from TinvTim, Urist & HookH"l999l) of the tele- 
scope+instrument system and averaged over the detector pixel 
size to obtain the observed surface brightness distribution. At 
this point, the derived model light profile can be directly com- 
pared with obse rved ones (for a deta iled description of the inver- 
sion and fit, see M arconi et al.ll2.Q03h . 



We considered the galaxy light profiles obtained by using 
the IRAF/STSDAS program ellipse to fit elliptical isophotes to 
the F160W continuum NICMOS images of NGC 3310, NGC 
4303 and NGC 4258. The ellipticity of the isophotes in the nu- 
clear region of all our galaxies (excluding the contribution from 
the central unresolved source, when present) is very close to 
zero (e = 0.15,0.05,0.0 for NGC 4303, NGC 3310 and NGC 
4258, respectively) allowing us to assume spherical symmetry 
(i.e. q = 1). In the case of NGC 3310, the isophote fit was limited 
to the region included by the circumnuclear ring visible in Fig. [2] 
in order not to be affected by the bright spiral arms present on 
large scale. The observed surface brightness profiles are shown 
in Fig. E 

We followed the above procedure to fit the observed surface 
brightness profiles within ~ 3" of the continuum peak. This is 
adequate to our goals since, as shown in the next section, the 
kinematical model fitting is performed only in the inner ~ 1". 
The best-fit parameters for all galaxies are presented in Tab. [3] 
where mass densities are computed for F = M/L = M /L 0i h, 
(its real value will be determined during the kinematical fitting 
procedure in Sec. [4j. The comparison between models and ob- 
served data is shown in Fig. [9] 

A nuclear point source was added for NGC 4303 and NGC 
4258 in order to account for the AGN emission. Indeed its pres- 
ence is clearly indicated by the little bump at log r 0.5, which 

represents the first Airy ring of the HST PSF. In the following 
analysis, we will assume that the central point sources are unre- 
solved emission from the AGN, and therefore they do not con- 
tribute to the total stellar light. An alternative possibility is that 
these central point sources are unresol ved star clusters common 
in late-type spirals dBoeker et al.l2004l) . However NGC 4303 and 
NGC 4258 are know n AGNs (see Sec[TJ. It has been shown by 
Ouillen et al. (2001) that unresolved infrared sources are found 
in the great majority of HST NICMOS images of Seyfert galax- 
ies and that their luminosities strongly correlate with both the 
hard X-ray and [OIII] line luminosity. This result strongly sug- 
gests a dominant AGN contribution to the IR emission. We thus 
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Fig. 9. Surface brightness profiles {erg cmT 1 s _1 arcsec^ 2 ) of the continuum emission from the F160W NICMOS images. The left, 
central and right panels refers to NGC 3310, NGC 4303 and NGC 4258, respectively. Dots with error bars (when visible) denote 
the observed data while the solid lines indicate the model fits (see text). The lower plots show the residuals after subtraction of the 
model fits. 



Table 3. Best-fit parameters for the stellar light density model. Mass densities (po, p5) are derived from luminosity density under 
the assumption that F = M Q /L h- X 2 re d 1S X 2 l^°f where dof represents the degrees of freedom of the fits (47 for NGC 4303, 37 for 
NGC 3310 and NGC 4258). 





q« 






a 


P 


y 


logPn" 


rf 


a* 


p 


1 PSF 


X ml 


NGC 4303 


1.0* 


2.28 


19.90 


2.35 


11.90 


2.00* 


-0.76 


8.51 


0.21 


7.87 


9.62 


13.09 


NGC 3310 


1.0* 


3.43 


0.66 


-0.71 


3.67 


2.00* 


3.40 


0.23 


1.79 


0.22 


4.40 


52.66 


NGC 4258 

a tw — a „ 


1.0* 


1.59 


4.55 

j i _ 


1.86 


0.14 


0.000* 












48.83 
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c In units of arcsec. 

d Point source flux in units of 10~ 17 erg cirT 2 s A -1 . 

* Parameter was held fixed, as explained in detail in 13.21 



consider unlikely that the point sources are due to unresolved 
star clusters. 

In NGC 4258, the observed brightness profile is well repro- 
duced by the spherical stellar light distribution of Eq. [2] and by 
a central point source. However in the cases of NGC 3310 and 
NGC 4303 a good match of the observed profiles required an 
additional component described by Eq. [2] but characterized by 
p* ,r* b ,a*,p*. 

Given p(r), one can easily derive M(r) the amount of stellar 
mass contained within the radius r. Fig.fTOlshows. for all galax- 
ies, the circular velocity due to the stars as a function of the dis- 
tance from the nucleus computed assuming M/L = Mq/L^h 
(solid line). The dashed line represents the circular velocity ex- 
pected from a BH with mass M BH = 10 7 M o . 

4. Kinematical Models 

The procedu re used to mode l the k inematical data was first 
described by Mac chetto et al.l ( 19971) and then refined by sev- 
eral authors dva n der Marel & van den Bosch II 19981: iBarfh et al] 
1200 It Marconi et al. 2003))- This procedure, described in detail 
in iMarconi et al.l (120031) . assumes that the ionized gas is circu- 
larly rotating in a thin disk placed in the principal plane of grav- 
itational galaxy potential. We neglect any hydrodynamical ef- 
fects like gas pressure. The gravitational potential is made by 
two components: the stellar one (whose mass distribution was 
determined in Sec. 13. 2\ and a dark mass concentration (the pu- 



tative BH), unresolved with the spatial resolution of our obser- 
vations. In computing the rotation curves we take into account 
the finite spatial resolution of observations, the intrinsic surface 
brightness distribution (which is the weight to average kinemat- 
ical quantities) and we integrate over the slit and pixel area. The 
free parameters characterizing the best fitting model are found 
by standard^ 2 minimization. 

In Sec. 14.11 we determine the Intrinsic Surface Brightness 
Distribution (ISBD h ereafter) which plays a very important role 
in the modeling (see IMarconi et al.ll2.Q05l for a detailed discus- 
sion). In Sec. [5] we describe the fitting of the rotation curves. 
IMarconi et al. (2005) have shown that the velocity dispersion 
depends on the adopted ISBD, therefore we decided to restrict 
the model fitting only to velocity. Moreover the observed veloc- 
ity dispersions of our galaxies are small (see Fig.|3]|4]and Fig. [8), 
not-mono tonic and likely affected by turbulent motions; we will 
discuss this issue further in Sec. [6] Summarizing, the free pa- 
rameters during the model kinematical fitting are: 

- the black hole mass Mbh\ 

- the stellar mass to light ratio M/L; 

- the coordinates of the BH position (X, Y) on the plane of the 
sky where the nuclear slit is centered at (0,0); 

- the galaxy systemic velocity V sys ; 

- the position angle of the disk line of nodes 0; 

The inclination i of the gaseous disk cannot be treated as a 
free parameter. Since the amplitude of the rotation curve scales 
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Fig. 10. Circular velocities expected from the stellar gravitational potential on the principal planes (dashed lines) with the as- 
sumption of unit M/L. The solid lines represents the circular velocities expected from the gravitational potential of a BH with 
M B h = 1O 7 M . The left, central and right panels refer to NGC 3310, 4303 and 4258, respectively. 



as v(r) 2 ~ M(r)sin 2 (i), i is tightly coupled with Mbh- The resid- 
ual, weak dependence of t he projected velocity on cos i (Eq. B2 
and B4 in Appendix B of iMarconi et al.ll2004l) can be used to 
derive i only when high S/N, int egral field kinematical data are 
available (e.g. Bart h et aTll2001l) but this is not the case of the 
data presented here. For this reasons we decided to work, for 
each galaxy, with an initial fixed value of the inclination and we 
then investigated the effects of varying i on the final best fit pa- 
rameters. 

Usually the PSF adopted in the modelling is the one at A 
6700 A directly obt ained with t he Tiny TIM m odeling software 
dKrist & Hookl 19991) . However. Ibressell (120061) pointed out that, 
by default, Tiny TIM computes the PSF with the presence of a 
Lyot stop which is not present when the G750M grating is used. 
The Lyot stop is mounted in front of the optical elements (e.g. 
G750L) close to the pupil plane and it supplements the focal 
plane mask used in coronographic mode by suppressing residual 
scattered and diffracted starlight from the collimated beam. It has 
a broadening effect on the PSF: with the Lyot-stop, in fact, the 
detector sees a smaller diameter primary mirror, has a broader 
Airy disk and a slightly more distant first Airy ring. A compar- 
ison of the radial profile of light for model imaging the PSF at 
6600 A with a nd without the Lyot stop, is shown in Fig. 1 of 
iDress el (2006). Therefore we followed the procedure suggested 
bv lDressell d2006l) to obtain from Tiny TIM the G750M PSF with- 
out the Lyot-stop. 

4.1. The Intrinsic Surface Brightness Distribution 

The observed kinematical quantities are averages over apertures 
defined by the slit width and the detector pixel size along the 
slit. The ISBD is thus a fundamental ingredient in the modeling 
of the kinematical quantities because it is the weight of the ave r- 
aging process (see, e.g., the discussion in Marc oni et alj |2005). 

iBarth et all (12001 1) derived the ISBD by deconvolving a con- 
tinuum subtracted emission line image and were successful in 
reproducing the microstructure of the rotation curves thus reduc- 
ing the overall^ 2 of their best fitting model. However, apart from 
the microstructure of the r otation curves which h as little impact 
on the final Mbh estimate. IMarconi et alj ([2003) showed that it 
is crucial to adopt a good ISBD in the central region where the 



velocity gradients are the largest. This has a strong impact on the 
quality of the fit and possibly on the final M BH estimate. When 
the observed line surface brightness distribution along the slit is 
strongly peaked in the nuclear region and has a profile which is 
a little different from that of an unresolved source, one cannot 
simply use a deconvolved emission line image but should try to 
reproduce the observed line surface brightness distribution with 
a parameterized intrinsic one and take into account the instru- 
mental e ffects. 

Recently IMarconi etail (120051) showed that the adopted ISBD, 
provided that it reproduces the observed one within the errors, 
does not affect the final BH mass estimate. In their paper, they 
showed that systematic errors on M B h due to the adopted ISBD 
were of the order of 0.05 - 0.08 in log Mbh- Notwithstanding the 
ISBD has an important effect on the quality of the velocity fit 
although not on the M B h value. 

We parameterized the ISBD with a combination of analytical 
functions. We used a combination of exponential, Gaussian and 
constant components (labeled respectively expo, gauss and const 
in tables of results) and their expression are, respectively: 

1(f) = r < r h 

I(r) = I oe - r/ro r> r h ( } 

j -d-n,) 2 

I(r) = ^2 e_ ^ W 
1(f) =0 r<r h 

l(r) = I r > r h {D > 

Io is the amplitude of the component and ro is its scale radius. 
The presence of a nuclear hole in the ISBD, modeled with a non- 
zero r/,, could be caused by varying ionization conditions of the 
gas (e.g. close to the nucleus the gas is too ionized to emit the 
line of interest) or by absence of the gas itself. In the previous 
equations other important parameters are 'hidden': these are q 
and 9 which are used to provide ISBDs with elliptical isophotes. 
q is the ratio between the minor and major axis of the isophotes 
and is the PA of the major axis. In practice: 

r = \x 2 + F 2 ] 1/2 (6) 
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where the coordinates (X, Y) are defined by: 
j X = (x - Xo) cos 9 + (y - y ) sin 6 

\Y = (-(*- xo) cos 6» + (y-yo) sin 0) 1 ; 

(*o,yo) represent the position of rotation center (the putative 
black hole) in a reference frame on the plane of the sky with x 
and y aligned along the North and East directions respectively. 

Summarizing the free parameters for each ISBD component 

are: 

- the amplitude Io; 

- the scale radius ro; 

- the hole radius r/, ; 

- the axial ratio q; 

- the position angle 6; 

- the coordinates of the symmetry center (xo,yo)', 

5. Kinematical results 

In this section we first determine the intrinsic surface brightness 
distribution for each galaxy and then we fit the observed velocity. 

5.1. Analysis of NGC 3310 

As explained in the Introduction, the galactic disk of NGC 3310 
is characterized by non-c ircular motions with str ong streaming 
along the spiral arms (e.g. Kregel & Sancisil200ll and references 
therein). Fig. [2] clearly shows that the spiral arms depart from a 
~ 6" ring, the starburst ring, and suggests that, beyond that ra- 
dius, motions might be significantly non-circular. Within the ring 
(which might indicate the corotation radius) the morphology of 
the emission line region is characterized by a peaked nuclear 
blob, about 1.2" in size, which appears unrelated with the star- 
burst ring. In order to avoid as much as possible obvious contam- 
ination by non circular motions we have decided to restrict our 
analysis to the nuclear blob, i.e. to the central 1.2"region. This 
region, which includes the high surface brightness peak shown 
in Fig.[3]is delimited by a dashed line. 

To determine the best ISBD we made several attempts 
with different compositions of the analytical functions described 
above. The best choices are the following: the one called A 
(composed of two Gaussian components, an exponential and a 
constant one) and the one called B (composed of two exponen- 
tial components, two Gaussian ones and a constant component). 
Their resulting best fitting parameters are shown in Tab. [4] and, 
in the left panel of Fig. Q~T] we compare the model with the ob- 
served flux distribution. 

As explained in Sec. [4] we fitted the velocities with fixed disk 
inclination and, as starting point for i, we adopted the inclination 
of the large scale galaxy disk (i = 40 deg). With both ISBD A 
and B, it is possible to reproduce the observed rotation curves 
(right panels of Fig.fTTI). Tab.[5]shows the corresponding best fit 
parameters. The reduced^ 2 of the fits are very large at variance 
with the apparently good agreement between data and models. 
The reason of this apparent contradiction is to be found in the 
small velocity errors, the formal ones of the line profile fitting. 
These errors are extremely small, of the order of 1 km/s, and cor- 
responding to 1/50 pixel. Although the S /N ratio of our spectra 
is large, such small errors values are not realistic and will be 
corrected as described below. 

The above kinematical fit was computed with the inclination 
i fixed at the value of the large scale disk inclination (z = 40 



deg ISanchez-Portal et al.ll2000l) . Recall that the choice of fix- 
ing the inclination is justified by its tight coupling with mass 
[v(r) 2 ~ M(r)sin 2 (i), see Sec. [4). We now verify the dependence 
of our model against the choice of this parameter. In particular 
we wish to check whether the quality of the fit might improve at 
different inclination. We thus fit the rotation curves with differ- 
ent fixed values of the inclination, from i = 5 deg to i = 85 deg, 
using the A flux distribution. Tab.[6]shows the best fit parameters 
obtained with the different inclinations. In order to establish the 
statistical significance of x 2 variations among different models 
we renormalized x 2 by adding in quadrature an extra velocity 
error to the data, Avo. The value of Avo is found by imposing 
that the model with the lowest^ 2 value has x 2 /d of = 1, where 
dof i s the number of degrees of freedom (see also lMarconi et al.1 
2004). With dof =30, the best fit models which are acceptable 
with a significance of 95% have X 2 scaled < 1-48. The depen- 
dence of x 2 on the adopted disk inclination is reported in the left 
panel of Fig. |T2j the filled square represents the inclination value 
i = 70 deg corresponding to the lowest^ 2 value. 

With the x 2 renormalization described above, all inclina- 
tions with i < 80 deg are acceptable at the 95% level and, ap- 
parently, our analysis produces an upper bound to the inclina- 
tion. However, the mass-to-light ratio increases sharply with de- 
creasing inclination and, indeed, also M/L scales with (sin;) -2 
as Mbh- We then considered the evolutionary synthesis models 
for stellar population derived by IVazquez & Leitherer I (2005) 
in order to derive a lower bound on i based on astrophysically 
acceptable M/L values. The typical stellar populations have 
H - K ~ 0.25 and the value of M/L monotonically increases 
with galaxy age reaching log M/L K < 0.6 i.e. log M/L H < 0.25 
(I Vazquez & Leitherer | |2005). Note that, for the sake of simplic- 
ity, we assume that the F160W band can be approximated with 
the standard H band. This physical upper limit on M/L can be 
translated into a lower limit on i which corresponds to an incli- 
nation of ~ 20 deg. To summarize, the constraints on inclination 
are that 20 deg < i < 80 deg. The model with the lowest^- 2 has 
i = 70 deg and the comparison between observed and model ro- 
tation curves is shown in Fig.[TTj Th e inclin ation value of i ~ 40 
deg, found by San chez-Portal et all d2000), refers to the large 
scale disk while the inner nuclear disk (which is clearly sepa- 
rated from the large scale structure, Fig. |2), is seen at a higher 
inclination. The best fit black hole mass for the 70 deg-model is 
log Mbh[M ] = 4.2 and, in general, all inclinations require simi- 
larly small black hole masses (log Mbh[Mq] < 5.6). Such values 
are unrealistically small and, indeed, we are only able to set an 
upper limit to Mbh- 

In order to estimate statistical errors associated with the best 
fit parameters, we present in the right panel of Fig. Q~2]the con- 
tours for the joint variation of Mbh and M/L. The grid has been 
computed with ISBD A and i = 70 deg, fixing the 2 parameters 
of interest (Mbh and M/L) and varying the others to minimize 
X 2 . Before computing confidence levels we have renormalized 
X 2 as described above so that the model with the lowest^- 2 value 
has x 2 reduced equal to 1 . Confiden ce levels for 2 parameters 
of interest are then found following lAvnil d 1976b : x 2 — x 2 • + 
2.3, 4.61, 6.17 for confidence levels 68.3%, 90% and 95.4%, re- 
spectively. Fig. Q~2] shows that it is not possible to measure the 
BH mass but only place an upper limit to it. By choosing a con- 
fidence level of 68.3%, the stellar mass to light ratio for this 
galaxy results M/L = O.47+° ^M /L 0jW which is the typical H- 
band value for spiral bulge s (0.6 + 0.2 in K, which correspond to 
0.3 ± 0.1 for H - K = 0.2: lMoriondo et al.lll998l) . For the black 
hole mass upper limit we chose a confidence level of 95% and 
the corresponding estimate results Mbh < 5.62 x 1O 6 M . 
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Table 4. NGC 3310: best fitting parameters of the two ISBDs adopted in the kinematical model fitting. 
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Fig. 11. NGC 3310. Left panel: best-fit models of the emission line surface brightness distribution compared with the observed data 
(dots with error bars) at the three different slit positions. The solid and dotted lines refer to ISBD A and B, respectively. From top 
to bottom the panels refer to NUC, OFF1 and OFF2 respectively. Central panel: best-fit models and observed data for the rotation 
curve obtained with A (solid line) and B (dotted line) flux distributions, assuming i = 40 deg. Right panel: best fit model of velocity 
obtained using ISBD A with the 'best' inclination (i - 70 deg) at the three different positions of the slit. 

Table 5. NGC3310: best fitting parameters of the kinematical model obtained with different ISBDs and i = 40°. 
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We then verified that at different inclinations the statistical 
errors on log M/L were roughly the same and that the upper lim- 
its on Mbh were scaling with (sin/) -2 as expected. By taking 
into account the allowed i variation (20 deg < i < 80 deg) the 
upper limit on the black hole mass varies in the range 5.0 x 10 6 
-4.2xl0 7 M o . 

In Fig. [13] we plot the expected velocity dispersion for the 
best fitting model at i = 70 deg which includes intrumental 
broadening and the contribution of unresolved rotation (dotted 
line). In order to match observations we added an intrinsic con- 
stant dispersion of 55 km/s (solid line). Overall data points do 
not show a deviation of more than 20 km/s from the model val- 
ues, after the intrinsic dispersion was taken into account. 

5.2. Analysis of NGC 4303 

The kinematics of NGC 4303 can be traced up to 4" from 
the galaxy n ucleus identified by the location of the near-IR 
point source dColina et al.l 12002). Although at 1 - 2" spatial 



resolution the motion within the central region (~ 10" of the 
nucleus) appears ordered and consistent with circular rotation 
dSchinnerer et al.1120021: IColina & Arribas| [l999) the presence of 
the inner bar with a siz e comparable with the spatial resolution 
(~ 2" size. lErwinl2004l) suggests caution in the kinematical anal- 
ysis. Indeed, the curve of the extended component shows irreg- 
ular behavior at about 3" from the nucleus (pix ~ 320 - 340 in 
Fig.|4|. Moreover, the Ha+[NII] surface brightness distribution 
along the slit clearly indicates the presence of a well defined nu- 
clear component as in the case of the other two galaxies in this 
paper and NGC 4041 dMarconi et alj|2003l) . This then suggests 
to limit our analysis to the nuclear component, i.e. to the central 
~ I" from the nucleus. 

We follow the same procedure adopted for NGC 3310. We 
considered different ISBDs but we found that only two of them 
are able to provide a reasonable fit of the observed data: the one 
called D, composed of three exponential components and a con- 
stant one, and the one called E, composed of three Gaussian 
components and a constant one. Fig. [14] shows the compari- 
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Table 6. NGC3310: effect of i variation on the best fit parameters and x 2 red - X 2 values have been rescaled by an extra error on 
velocity Av . All models adopt ISBD A. 
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Fig. 12. NGC 3310. Left: dependence of^ 2 on the adopted disk inclination. Horizontal lines indicate the 95 %, 90% and 68% 
confidence levels. The filled square represents the point corresponding to a^ 2 escaM = 1 for ;' = 70 deg. The^ 2 value at ;' = 85 
deg is out of the plot range. Right: x 2 contours for the joint variation of Mbh and M/L (with i = 70 deg). Contour levels are for 
X 2 - X 2 m i n + 2.3,4.61,6.17 corresponding to 68.3%, 90% and 95.4% confidence levels. The dots indicate the M BH and M/L values 
for which the x 2 minimization was actually computed. 



Table 7. NGC 4303: best fitting parameters of the two ISBDs adopted in the kinematical model fitting. 
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Fig. 13. NGC 3310. Expected velocity dispersions for the best fitting rotating disk model with i = 70 deg compared with observations 
(dots with error bars) assuming a constant intrinsic dispersion of 55 km/s (solid line) compared with the dispersion expected from 
pure rotation (dotted line). 
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Fig. 14. NGC 4303. Left panel: best-fit models of the emission line surface brightness distribution compared with the observed data 
(dots with error bars) at the three different slit positions. The solid and dotted lines refer to ISBD D and E, respectively. From top 
to bottom the panels refer to NUC, OFF1 and OFF2 respectively .Center panel: best-fit models (solid and dotted lines) and observed 
data (dots with errorbars) for the rotation curves at the three slit positions. The solid and dotted lines denote models obtained with 
the D and E ISBD, respectively. The models were computed assuming i = 20 deg. Right panel: as in the left panel, but the solid line 
is the model computed for the best inclination i = 70 deg and ISBD D. 



son between observed and model surface brightness distributions 
while the best fit ISBD parameters are placed in Tab. [7] 

As for the case of NGC 3310, we started the velocity fitting 
with a fixed inclination value, i = 20 deg corresponding to the 
large scale value one, and we used both ISBDs described above. 
The nuclear rotation curves of NGC 4303 have a less regular 
appearance than NGC 3310 and the agreement with the mod- 
els might thus be worse. A blind x 2 minimization can result in 
non-physical parameters values. In particular the position of the 
kinematical center of rotation can move away by > 0.5" from the 
continuum peak which clearly identifies th e galaxy nucleus (e.g. 
ISchinnerer et alj|2002h IColina et alj|2002l) . To avoid this prob- 
lem we have constrained the x,y coordinates of the kinematical 



center to be located within +0.2" of the continuum peak posi- 
tion. The best fits of the rotation curves with ISBDs D and E are 
shown in the left panel of Fig. [14]and the best fitting kinematical 
parameters are shown in Tab. [8] 

ISBD D provides a much better fit of the observed velocities 
and the disagreement in the best fit parameter values among the 
two models is readily explained by the difference in reduced x 2 
values. Model D has two unsettling characteristics: an unphysi- 
cally high M/L (which might be due to the low disk inclination) 
and a PA of the line of nodes which is very close to 90 deg, a ge- 
ometry in which only a small fraction of the rotational velocity 
is projected along the line of sight. 
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Table 8. NGC 4303: best fitting parameters of the kinematical model obtained with different ISBDs and i = 20°. 
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Table 9. NGC 4303: effect of i variation on the best fit parameters and x 2 ed - X 2 values have been rescaled by an extra error on 
velocity Av . All models adopt ISBD D. 
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Again, we have investigated whether different inclinations 
can provide a better fit to the observed rotation curves. Tab. [9] 
shows the resulting best fit values obtained for different incli- 
nations from i = 5 deg to i = 85 deg, using the D ISBD. As 
previously, x 2 have been rescaled so that the model with the 
lowest x 2 has a value of 1 for the reduced x 2 - The x 2 depen- 
dence on the adopted disk inclination is shown in the left panel 
of Fig. [15] where the model with the lowest x 2 is marked with 
a filled square. In the case of NGC 4303, 24 degrees of freedom 
imply that the 95% confidence level is found for^? eJcoferf = 1.52. 
Above that value all models are to be discarded with a 95% sig- 
nificance implying i > 40 deg. The model with the lowest x 2 
has i = 70 deg, M BH = 4.90 x 10 6 M o and is plotted in the right 
panel of Fig. [14] The first modeling attempt, made at i = 20 deg, 
should be discarded and the non-physical value of M/L is not 
problematic. Moreover, at high inclination values the PA of the 
line of nodes is sign ificantly different from 9 de g and is similar 
to the one found by Colina & Arribasl d 19991) and lKoda &~S ofue 
(2006) for the arcsec scale gas rotation (130 - 180 = -50 vs 130 
and -44.5, respectively). As in the case of NGC 3310, the incli- 
nation of the nuclear disk is smaller than that of the large scale 
galactic disk (z = 45 deg). 

Apart from the consistent kinematical parameters found, the 
agreement between model and observed rotation curve is good 
except for OFF1 where discrepancies might indicate the pres- 
ence of non-circular motions. 

In order to estimate statistical errors associated with the best 
fit parameters, we present in the right panel of Fig. [15] the con- 
tours for the joint variation of Mbh and M/L. The grid has 



been computed with ISBD D and i = 70 deg, fixing 2 interest- 
ing parameters (Mbh and M/L) and varying the others to min- 
imize^ 2 (see Sec. 15.11 ). By choosing a 68.3% confidence level, 
M BH = (5.0)!^ x 10 6 M o and M/L = (O.1)!^M / L °.", in- 
dicative of a veryyoung stellar population in the central region 
(IColina et alj|2002b . 

Fig.fTBIprovides an alternative visualization of the best fitting 
model. The velocity field weighted with the intrinsic line flux 
distribution is shown prior to the averaging over slit and pixel 
sizes. The velocity asymmetries with respect to the lines of nodes 
(dashed line) are due to the intrinsic line flux distribution, which 
is not symmetric with respect to the BH position. 

To check for the robustness of the BH detection we tried to 
fit the rotation curves without a BH but only with the gravita- 
tional potential due to the stars. We worked at i = 70 deg and 
with D ISBD. The result of this fit is shown in the left panel of 
Fig- El At first sight the no-BH fit does not appear significantly 
worse than the one with a BH but it is nonetheless statistically 
unacceptable. After rescaling velocity errors with 6v = 6.9km/ s, 
as described above, the 68.3%, 90%, 95.4% confidence levels are 
X 2 = 27.8, 30.1 and 31.7, respectively while the model without 
the BH has x 2 = 33.8. However, given the non perfect agree- 
ment between our model and the data, we do not consider our 
Mbh detection completely reliable. 

As in the case of NGC 3310, we verified that at different 
inclinations the statistical errors on log M/L and log Mbh were 
roughly the same. The allowed range of inclinations (i > 40 deg) 
is then translated into a possible range for Mbh which can vary 
between 6.0 x 10 5 and 1 .6 x 1O 7 M where the extreme have been 
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Fig. 15. NGC 4303. Left: dependence of x on the adopted disk inclination. Horizontal lines indicate the 95 %, 90% and 68% 
confidence levels. The filled square represents the point corresponding to ax 2 , , = 1 for i - 70 deg. Right: x 2 contours for the 
joint variation of M BH and M/L (with i = 70 deg). Contour levels are for^- 2 = x 2 ■ + 2.3,4.61, 6.17 corresponding to 68.3%, 90% 
and 95.4% confidence levels. The dots indicate the Mbh and M/L values for which the^ 2 minimization was actually computed. 
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Fig. 16. Map of the velocity field in the best fitting model of NGC 4303. The velocity field is weighted with the intrinsic line flux 
distribution and is shown prior to averaging over slit and pixel sizes. The dashed areas denote the three parallel slits. The filled 
square represents the position of the black hole and the dashed line is the line of nodes. 



computed by taking into account also the 95% confidence errors 
on \ogM BH (-0.4, +0.15). 

In Fig. [T7J we plot the expected velocity dispersion for the 
best fitting model at i = 70 deg which includes intrumental 
broadening and the contribution of unresolved rotation (dotted 
line). In order to match observations we added an intrinsic con- 
stant dispersion of 50 km/s (solid line). Overall data points do 



not show a deviation of more than 20 km/s from the model val- 
ues, after the intrinsic dispersion was taken into account. 

5.3. Analysis of NGC 4258 

An inspection of Fig. [8] shows that the available kinematical pa- 
rameters for this galaxy are limited to the nuclear component i.e. 
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Fig. 17. NGC 4303. Left panel: best-fit model without a BH at i = 70 deg and ISBD D. The dotted line represents the best fit model 
with a BH, obtained with the same inclination and ISBD. Notation as in the Fig. [14] Right panel: expected velocity dispersions for 
the best fitting rotating disk model with i = 70 deg compared with observations (dots with error bars), assuming a constant intrinsic 
dispersion of 50 km/s (solid lines). The dotted line is the velocity dispersion expected only by unresolved rotation. 



within 0.5" of the continuum peak. We followed the same proce- 
dure outlined for NGC 3310 and NGC 4303 with the difference 
that, in this case, we fitted simultaneously ours and archival data. 
Effectively we fit 6 slit positions at the same time. 

The starting assumption is that the nuclear slits are centered 
on the same point, as confirmed from the fact that the broad Ha 
components have very similar fluxes in the two nuclear spectra 
(see Fig. [6}. 

We first tried to fit the observed surface brightness distri- 
bution with different combinations of analytical functions, but 
we could find only one ISBD that provided a good agreement 
with the data. Left panel of Fig. [18] and Tab. [10] show the best 
fit model and the parameter values, respectively. It is clear from 
the figure that the surface brightness drops by a factor ~ 10 from 
the nuclear to the offset positions. Note that the observed surface 
brightness in the NUC slit of the archival spectra is a factor 2 
larger than in our data. This is simply due to the fact that our 
data use the 0.2" slit instead of the 0.1" one. 

To fit the observed rotation curves, we started assuming that 
disk inclination is ; = 80° as in the much smaller maser disk 
dHerrnstein et alj[l999h . The best-fit model is shown in the right 
panel of Fig. [18] and the best fit parameter values are displayed 
in Tab. E] 

As in the cases of NGC 3310 and NGC 4303, we investigate 
the influence of the disk inclination on quality of the fit. Tab. [12] 
shows the best-fit parameters obtained at different values of the 
inclination. 

The model with the lowest x 1 value is the one at i = 60 deg 
and, after rescaling, all the models have reduced y 1 close to 1 . 
The left panel of Fig. [19] shows the dependence of x 2 on the 
adopted disk inclination (the filled square marks the model with 
the lowest^ 2 ). In this figure we did not draw the 90% and 95% 
confidence level because they are larger than all the values ofx 2 
found by our fitting. This means that its is not possible to place 
any constraints on the inclination value but it is interesting to 
note that the best inclination value is remarkably si milar to the 
inclin ation of the galactic disk which is =s 64 deg dCecil et alj 
fl992h . 



Finally, we present in the right panel of Fig.[l9]the contours 
for the joint variation of Mbh and M/L computed for the disk 
inclination of i = 60 deg. The contour shapes allow to estimate 
M B h, but only an upper limit for the M/L ratio. At the 68.3% 
confidence level, the black hole mass is M BH = (7.9)^x 10 7 M o 
while the 95% upper limit on M/L is M/L < 1 .SMq/Lqjj. 

In the left panel of Fig. [20]we compare the observed kinemat- 
ics with the models with (solid line) and without (dashed line) a 
BH computed at i = 60 deg. As in the case of NGC 4303, the no- 
BH fit does not appear much worse than the one with a BH but 
is statistically unacceptable. After rescaling velocity errors with 
6v = 32.9km/ s, as described above, the 68.3%, 90%, 95.4% con- 
fidence levels are x 2 = 35.3, 37.6 and 39.1, respectively while 
the model without the BH has x 1 = 57.8. Moreover, the no-BH 
model is characterized by a very high value of the mass-to-light 
ratio M/L = 5.O8M /L Oj h which is larger than the maximum 
predicted by normal stellar populations (see the discussion in 
the case of NGC 3310 in Sec. 5.1) indicating the presence of a 
dark mass. Combining these two facts we consider that with our 
data we have reliably detected the presence of the BH in NGC 
4258. 

Our black ho le mass value i s com patible with the H2O- 
maser estimate by Gree nhill et al.l (Il995allbl) and lHerrnstein et ail 
dl999l) which is M BH = (3.9 + 0.1)x 10 7 M o at slightly above the 
lcr level. Given the poor match between the data and a simple 
pure rotating disk model (i.e. no streaming or warps), this must 
be nonetheless considered an encouraging result. 

As in the previous cases, we verified that the statistical er- 
rors on log M/L and log Mbh do not vary significantly with the 
adopted inclination. However, in this case we could not place 
any constraint on the inclination thus, in principle, we could only 
provide a lower limit on the BH mass value. The nuclear gas disk 
we are mapping is spatially intermediate between the large scale 
galactic disk and the maser disk. The two disks have different in- 
clinations ~ 60 and 80 deg, respectively and it is thus reasonable 
to assume that the inclination of the nuclear gas disk, connecting 
the large to the small galactic scales, is limited in the 60 - 80 deg 
range. With this assumption the black hole mass from gas kine- 
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Fig. 18. NGC 4258. Left panel: best-fit model (solid line) of the emission line surface brightness distribution compared with the 
observed data (dots with error bars) at the different slit positions. From top to bottom the panels refer to NUC, OFFl and OFF2 
respectively. The different columns denote our and archival data. Right panel: best-fit model (solid line) and observed data (dots 
with error bars) for the rotation curves. The model is computed for i = 80 deg. The meaning of the columns is the same as in the 
left panel. 

Table 10. NGC 4258: best fitting parameters of the ISBD adopted in the kinematical model fitting simultaneously our and archive 
data. 



Id function i 



1" 



roiC) nA") Xot(") y 0i (") 



expo 
expo 
const 



14139 
1173 
2 



0.01 
0.03 



0.00 
0.00 



0.04 
0.07 



-0.01 
-0.04 



-0.1 
-0.4 



0.49 
0.30 



" Units of erg s 

Table 11. NGC 4258: best fitting parameters from the velocity fit for i = 80 deg fitting simultaneously our and archive data. 



Flux" x (") yo(") logM'^ log Mj L c 6Q Kn,(km/s) X L,~ 



0.02 -0.01 



8.02 



0.70 



235.44 



466.45 



34.96 



" Adopted flux distribution 

b Units of M 

c Units of M Q /L BM 



matics is confined in the range 2.5 x 10 7 - 2.6 x 1O 8 M where 
the extreme have been computed by taking into account also the 
95% confidence errors on log M BH (-0.5, +0.4). 

In the right panel of Fig.|20]we plot the expected velocity dis- 
persion for the best fitting model at i = 60 deg which includes 
intrumental broadening and the contribution of unresolved ro- 
tation (solid line). In contrast with NGC 3310 and NGC 4303, 
the observed velocity dispersion can be well matched with unre- 
solved rotation without the inclusion of any intrinsic dispersion. 

The non-perfect agreement between our model and the data 
suggests that rotation curves might be somehow affected by non 
circular motions. Indeed a jet-cloud interaction is responsible for 
the hot X-ray coco ons which give rise to the anomalous arms 
jWilson et al.ll200il) . Their position angle, i.e. the region where 



the X-ray cocoons are interacting with the dense interstellar gas, 
is ~ -33 deg, very close to the slit PA of the archival data 
(~ 146 = -34). Since our data are at an almost perpendicu- 
lar slit PA and likely less affected by non-circular motions, we 
have tried to perform the fit without the archival data. We found 
no significant change in the black hole mass estimate although 
a slightly better agreement between model and observed veloc- 
ities. This indicates that the non-circular motions do not affect 
the final black hole mass estimate. 

6. Discussion 

In Sec. 15.1115. 21 and l53l we have presented the kinematical mod- 
eling of the three Sbc spiral galaxies analyzed in this paper, NGC 
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Table 12. NGC4258: effect of i variation on the best fit parameters and x 1 a- X 1 values have been rescaled by an extra error on 
velocity Avq. 



K°) 


\ogM a BH log M/L" 


V. m (km/s) 


on 


(y 2 Y 

^ rescaled ' 






Fit of velocity (Av = 33.74/0?™"')" 




5 


9.69 


-0.74 


465.60 


239.22 


l.ll 


10 


9.12 


-1.21 


466.31 


239.84 


1.09 


15 


8.78 


-1.15 


466.78 


239.54 


1.09 


20 


8.56 


-2.10 


461.66 


238.51 


1.10 


25 


8.37 


-1.80 


465.01 


238.95 


1.08 


30 


8.26 


-1.92 


464.67 


238.98 


1.07 


35 


8.15 


-1.74 


466.36 


238.68 


1.06 


40 


8.06 


-1.83 


465.47 


238.35 


1.05 


45 


8.03 


-2.24 


463.81 


237.64 


1.04 


50 


7.99 


-2.22 


465.32 


237.07 


1.03 


55 


7.97 


-1.96 


465.13 


236.85 


1.02 


60 


7.96 


-2.30 


467.98 


237.04 


1.00 


65 


8.00 


-1.70 


468.68 


236.14 


1.00 


70 


8.01 


-1.64 


472.64 


236.69 


1.00 


75 


8.09 


-1.76 


472.40 


236.45 


1.02 


85 


8.47 


-2.07 


497.09 


244.23 


1.12 



8 Units of M . 

b Units of M e /L QM . 

c Rescaled^ 2 with errors computed as Av ( 2 =Av 2 + Av 2 ,. 
d Systematic error adopted to renormalize^ 2 . 




Fig. 19. Left panel: Dependence of y 1 on the adopted disk inclination during NGC4258 fitting only our data. Dotted line indicates 
the 68.3% confidence levels (lines indicating respectively 90% and 95.4% would have value 1.45 and 1. 34). Right panel: For the 
same galaxy,^ 2 contours for joint variation of Mbh and M/L using the best ISBD and i = 60 deg. Contours levels are for^ 2 = 
Xmin + 2.3, 4.61, 6.17 corresponding to 68.3%,90% and 95.4% confidence levels. 



3310, NGC 4303 and NGC 4258. We first worked with a fixed 
value of inclination and then we found the range of possible val- 
ues for M BH taking into account the allowed values of inclina- 
tions. Previous works on spiral galaxies have made the simplified 
assumption that nuclear gas disks used for Mbh detection have 
the same inclination as the large scale galactic disks. In fact, in 
contrast to elliptical galaxies, it is not possibile to determine the 
gas disk morphologies and inclinations from HST images of spi- 
ral galaxies. In this paper we have considered the effects of the 
unknown inclination which results in a range of Mbh values or 
upper limits. The M B h values derived for a given disk inclination 
are characterized by statistical errors as shown in the previous 
section but the same Mbh values vary as (sin iy 2 . Thus a range 
of allowed inclinations results in a range of possible BH masses. 



The observed rotation curves of NGC 3310 within the in- 
ner 0.6" from the nucleus (r < 50 pc) agree very well with 
a circularly rotating disk model. The observed velocity disper- 
sions cannot be matched with a simple rotating disk but re- 
quire a constant intrinsic velocity dispersion of ~ 55 km/s. This 
might indicate the presence of turbulence, non-circular or non- 
gravitational motions. The former would just contribute to the 
dispersion of the gas but not to the bulk motions and therefore 
would not invalidate the gas kinematical method. Conversely, 
the latter ones could invalidate the measurements depending on 
their strength co mpared to the ordered circular rotation (see, e.g. , 
the discussion in Barth et al. 2001; Verdoes Kleiin et al. 2006; 
iMarconi et al.ll2.004l) . However, the good agreement between the 
observed and model rotation curves suggests that the contribu- 
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Fig. 20. NGC 4258:Le/f panel: best-fit model (solid line) computed for i = 60° fitting simultaneously our and archive data (dots 
with error bars). From top to bottom the panels refer to NUC, OFF1 and OFF2 respectively. The dashed line is the best fit model 
obtained without a BH; even if the difference between the two cases does not appear significative, the model without a BH provides 
an unphysically high value for stellar mass to light ratio which is M/L = 5.O8M /L 0i h- Right panel: expected velocity dispersions 
for the best fitting rotating disk model with i = 60 deg (solid line) compared with the observations (dots with error bars). 



tion of non-circular or non-gravitational motions is small and 
does not invalidate the measurement. From the analysis of the 
rotation curves it is not possible to find evidence for the exis- 
tence of any BH since they can be accounted for by the stel- 
lar mass. The derived mass to light ratio (for the best model 
at i = 70 deg) is M/L = (0.47) + gffM Q /L Q , H and is consis- 
tent w ith the typical numbers for spiral bulges (Mor iondo et alj 
Il998l) confirming the validity of the analysis. The spatial resolu- 
tion of our observations is about 0.07"(FWHM of the STIS PSF 
at 6700A) and within a sphere with that diameter the stellar mass 
is M» ~ 1O 7 M . This fact, combined with the upper limit on the 
BH mass (M BH < 4.2 x 1O 6 M , at the 96% level for the lowest 
inclination allowed), indicates that we were not able to detect the 
BH because our spatial resolution is not sufficiently high. 

The situation for NGC 4303 is different with respect to NGC 
3310. The rotation curves require the presence of a BH with 
mass M BH = 5.0+°^ x 1O 6 M (6.0 x 10 5 - 1.6 x 1O 7 M when 
taking into account i > 40 deg) but the agreement between ob- 
served and model velocities is not as good especially at OFF1. 
With the inclusion of an intrinsic velocity dispersion of 50 km/s, 
the observed velocity dispersions in the offnuclear slits are in fair 
agreement with the model ones while in the nucleus they show 
deviations of the order of 10 km/s. Although the model without a 
black hole is statistically worse than that with a BH, the combi- 
nation of the above facts does not allow us to consider the black 
hole mass estimate in NGC 4303 completely reliable. 

For NGC 4258, the fair agreement between observed and 
model rotation curves might indicate the presence of non- 
circular motions but the observed velocity dispersion is com- 
pletely accounted for by unresolved rotation. These two facts 
alone would not allow us to say anything about the reliability 



of the BH mass estimate. However, NGC 4258 was chosen as 
a benchmark for the gas kinematical method because it rep- 
resents the second best case for a BH (after the Milky Way), 
with an accurate mass e stimate (Mbh - (3.9 + 0.1) x 1O 7 M ; 
iGreenhill et al.l 1995a b; Herrnstei n et al. The agreement 

between our mass estimate 1O 7 M ) (z = 60 deg) and the 

//2<9-maser value (the difference is just slightly above lcr) indi- 
cates that the gas kinematical method is reliable even if the ob- 
served velocities present deviations from the circularly rotating 
disk model. We now compare the BH mass estimates (or upper 
limits) in spiral galaxies with the known correlations between 
Mbh and host galaxy structural parameters. Since our work has 
doubled the number of gas kinematical measurements in late 
type spiral galaxies it is worthwhile to consider the sample of 
all late type spirals studied so far. 

This is shown in Tab. [T3l where for all galaxies studied in our 
project and those available from the literature we report the mea- 
sured BH masses or upper li mits and the ratio with th e expected 
values from the M BH - L S Dn dMarconi & Hunll2003t eq. [8]) and 

M BH I a, dTremaine et alJl2002t iFerrarese & Fordl 120051; eq. 
andfTUb correlations: 

log(AWM ) = 8.21 + 1.13(log(L^/L 0jX ) - 10.9) (8) 

log(M BH /M Q ) = 8.13 + 4.02 log(o-,/200 km/s) (9) 

log(M BH /M ) = 8.22 + 4.861og(cr c /200£m/s) (10) 

BH masses are taken from thi s work (NGC 3310, NGC 
4303 ) from lAtkinson et alj 12001 (NGC 1300, NGC 2748), 
from iMarconi et alj 120031 (NGC 4041) and from the compi- 
lation by Ferr arese & Fordl 120051 (Milky Way, NGC 4258 and 
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Table 13. Comparison between Mbh estimates in late type spirals with the expectations from the BH-spheroid scaling relations. 
The intrinsic scatter of the correlations is ~ 0.3 in log Mbh and the uncertainty on the determination of the luminosity is 0.3 dex 
(see text). 



Galaxy 


Type 


Tf 


l°g AC 




log M B „/M BH (L K y 


CT C ' 


log M BH /M BH (a c y 




log M BH IM BH (cr e y 


NGC1300 


SB(rs)bc 


18.8 


1 8+0-2 
'• 8 -0.2 


10.0 


0.6 


90 


1.3 


87 


1.1 


NGC2748 


SAbc 


23.2 


1 6 + °' 2 


9.8 


0.6 


79 


1.4 


83 


1.0 


NGC3310 


SAB(r)bc 


17.4 


<7.6 


9.6 


<0.9 


101 


< 0.8 


84 


< 1.0 


NGC4041 


SA(rs)bc 


19.5 


<7.3 


9.7 


< 0.4 


92 


< 0.7 


88 


<0.6 


NGC4303 


SAB(rs)bc 


16.1 


6-6^ 


10.2 


-0.8 


108 


-0.3 


84 


0.0 



NGC4258 


SAB(s)bc 


7.2 


7 cq+0.04 
'• J -0.04 


10.3 


0.09 


120 


0.45 


148 


-0.19 


Milkyway 


Sbl-II 


0.008 


6 60 +0 03 


10.2 


-0.9 


100 


-0.2 


100 


-0.3 


M81 


SA(s)ab 


3.9 


7 84 +0 1 

'•°^-0.07 


11.0 


-0.4 


174 


-0.1 


165 


0.0 



" units of Mpc. 
b units of M Q . 

c units of L Q K . 

d M BH (L K ) is the value expected from the Marco m~& Huntl |2003) correlation. 

' units of km/s. 

f M B n(cr c ) is the value expected from the Ferrarese & Ford (2005) correlation. 
s M BH (cr e ) is the value expected from the Tremai ne et alJ d2002h correlation. 



M81) to which we refer for the proper references for each 
galaxy. K band l umino sities of the spheroid are taken from 



Marconi & Huntl 0003 (M ilky Way, NGC 4258, M81), from 
Dong & De Robertisl (l2006h (NGC 2748 and NGC 4041). For 
the other galaxies (NGC 1300, NGC 3310 and NGC 4303) we 
took the total K b and magnitude fr om the 2MASS extended 
source catalogue dJarrett et alj 2000[) and we appl i ed the near- 
IR bulge-total correction by iDong & De Robertisl d2006l) . The 
correcting factor depends on the galaxy morphological type T 
and for our galaxies (T = 4.0), it is Am = M sp h - M lo , = 2.0. 
We then transformed the relative to absolute magnitudes using 
the distances shown in the table and we obtained the luminosi- 
ties using the solar K band magnitude M^ = 3.28. The uncer- 
tainty on the determination of /^-band luminosity of each galaxy 
is 0.3 dex and is due to the observed dispersi on on the bulge-total 
corre ction for a given morphological type dDong & De Robertisl 
2006). The stellar velocity dispersions required by the corre- 
lations in Eq. [9] and [10] are different: Trema ine et alj d2002l) 
use the luminosity weighted velocity dispersion cr e measured 
within r e i.e. with i n the bulge half-light radius. On the contrary, 
Ferrarese & Ford (2005) use the central velocity dispersion cr c , 
i.e. measured within an aperture of r — |r c . For all galaxies ex- 
cept the Milky Way an d M81 we used the values measured by 
Batcheldor et a Ud2005b . O therwise we used t he compilations by 
Tremaine et aiTd2002l) and lFerrarese & Fordl d2005l) . 



When inspecting the logarithms of the ratios between the 
measured BH masses and the expectations from the correla- 
tions one has to bear in mind that the intrinsic dispersion of the 
correlations is ~ 0.3 in log Mbh (M arconi & Hunjl2003T) . The 
upper limit on the BH mass in NGC 3310 is consistent with 
all correlations since it is larger than all expected Mbh values. 
The BH mass measurement in NGC 4303 is less reliable, as 
discussed above. However, it is comforting that the measured 
Mbh value is in very good agreement with the Mbh - cr corre- 



lations, if their intrinsic dispersion is taken into account. Mbh 
is a factor of 6 smaller than expected from the Mbh - Lk cor- 
relation however it should be recalled that the bulge luminosity 
of NGC 4303 has been me asured using the scaling relation by 
Dong & De Robertii d2006l) and the observed dispersion for Sbc 
spiral galaxies is of the order of 0.3 in log Lk- 

In Fig|2T] we compare spiral galaxies (filled squares) with 
the known Mbh i '^so IixMbhi 'o 'e and Mbh/o 'c correlations from 
the co m pilations by iMarconi & Hunl d2003l) : iFerrarese" & Fordl 
(2005J); iTremaine et al.1 d2002l) . The empty squares represent 
early type galaxies and, for simplicity, we have not shown er- 
ror bars. 

Considering that each correlation has an intrinsic dispersion 
of 0.3 in log Mbh we can conclude that spiral galaxies do not 
show evident deviations from the correlations followed by early 
type galaxies, apart for a slightly larger scatter. A statistical anal- 
ysis on the behavior of spiral galaxies has to wait for a larger 
number of BH mass detections. 

In conclusion, in order to show the difficulty of this kind of 
measurements, we compute the diameters of the BH spheres of 
influence (2 fg ) for each galaxy and compare them with the 
available spatial resolution (FWHM ~ 0.07 arcsec). These are 
< 0.08, 0.08 and 0.44" for NGC 3310, NGC 4303 and NGC 
4258, respectively. Except for NGC 4258, which is a nearby spi- 
ral galaxy, the black hole spheres of influence are at the limit 
of the HST spatial resolution, indicating how challenging these 
measurements are. 



7. Conclusions 

We have presented new HST Space Telescope Imaging 
Spectrograph (STIS) observations of three Sbc galaxies: NGC 
3310, NGC 4303 and NGC 4258. We used Ha A 6564A, [Nil] 
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Fig. 21. Correlations between Mbh and the host galaxy luminosity (right panel), the central stellar velocity dispersion, cr c (middle 
panel) and the effective velocity dispersion cr e (left panel). The empty squares denote early type galaxies (E/SO) from the literature 
while the filled ones refer to the spiral galaxies in TabfT3l 



AA 6549, 6585 A and [SII] AA 6718, 6732A to study the kinemat- 
ics of nuclear ionized gas disks at ~ 0.07" spatial resolution. 
Our NGC 4258 data were analyzed in conjunction with available 
HST/STIS archival data. In all galaxies, we derived the kinemat- 
ical moments by fitting a Gauss-Hermite expansion to the ob- 
served line profiles and, in the case of NGC 4258, we also took 
into account the presence of a broad Ha component. 

We have modeled the observed kinematics by assuming that 
the gas is confined in a thin circularly rotating disk where hydro- 
dynamical effects can be neglected. The disk rotates in a gravi- 
tational potential which is due to the galaxy stars and to a black 
hole. For each galaxy the stellar mass distribution was derived 
from HST-NICMOS images by fitting the observed stellar sur- 
face brightness. In computing the model kinematics, we have 
taken into account the finite spatial resolution of the observations 
and we have averaged over the slit and pixel area, weighting with 
the intrinsic surface brightness distribution of the emission lines. 
The free fit parameters were determined with a^ 2 minimization. 

In NGC 3310, the observed rotation curves are well repro- 
duced with the circularly rotating disk model. However, the spa- 
tial resolution of the observations only allowed us to set an upper 
limit to the BH mass, M BH < 5.6 x 1O 6 M , for a disk inclina- 
tion of ; = 70 deg. When taking into account the allowed disk 
inclinations, M BH varies in the range 5.0 x 10 6 - 4.2 x 10 7 M o . 

In NGC 4303 the rotation curves require the presence of a 
black hole with mass M BH = (5.0)+^ x 1O 6 M (z = 70 deg) 
but the agreement with the circularly rotating disk model does 
not allow to consider this detection completely reliable. If the 
allowed inclination values are taken into account, M B h varies in 
the range 6.0 x 10 5 - 1.6 x 10 7 M o . 

Finally, in NGC 4258, the simultaneous modeling of our and 
archival data constrains the BH mass to M BH = (7.9+35) x 10 7 
(i = 60 deg) and, taking into account reasonable limits for the 
inclination, M BH is in the range 2.5 x 10 7 - 2.6 X 10 S M G . The ob- 
served rotation curves are not well matched by the model but the 
good agreement wit h the accurate i^O-mas er estimate (M B h - 
3.9± 0.1) x 10 7 MM lGreenhill et all fl995allbT) : lHerrnstein et alj 
1999)) is an important indication of the reliability of the gas 
kinematical method. 



We have combined the BH mass estimates in late type spirals 
presented in this paper with others from the literature. Although 
the sample is still small (6 detections and 2 upper limits), spiral 
galaxies seem to follow the same BH-spheroid scaling relations 
as early type galaxies. 

The comparison between the BH spheres of influence and the 
spatial resolution of HST observations underlines the difficulties 
connected to the study of BH in late type spiral galaxies and in- 
dicates that higher spatial resolution is required for a significant 
step forward. 
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